Bayesian methods of astronomical source extraction

We present two new source extraction methods, based on Bayesian model selection and using the Bayesian Information Criterion (BIC). The first is a source detection filter, able to simultaneously detect point sources and estimate the image background. The second is an advanced photometry technique, which measures the flux, position (to sub-pixel accuracy), local background and point spread function. We apply the source detection filter to simulated Herschel-SPIRE data and show the filter's ability to both detect point sources and also simultaneously estimate the image background. We use the photometry method to analyse a simple simulated image containing a source of unknown flux, position and point spread function; we not only accurately measure these parameters, but also determine their uncertainties (using Markov-Chain Monte Carlo sampling). The method also characterises the nature of the source (distinguishing between a point source and extended source). We demonstrate the effect of including additional prior knowledge. Prior knowledge of the point spread function increase the precision of the flux measurement, while prior knowledge of the background has onlya small impact. In the presence of higher noise levels, we show that prior positional knowledge (such as might arise from a strong detection in another waveband) allows us to accurately measure the source flux even when the source is too faint to be detected directly. These methods are incorporated in SUSSEXtractor, the source extraction pipeline for the forthcoming Akari FIS far-infrared all-sky survey. They are also implemented in a stand-alone, beta-version public tool that can be obtained at http://astronomy.sussex.ac.uk/$\sim$rss23/sourceMiner\_v0.1.2.0.tar.gzComment: Accepted for publication by ApJ (this version compiled used emulateapj.cls

Identification of new transitional disk candidates in Lupus with Herschel

New data from the Herschel Space Observatory are broadening our understanding of the physics and evolution of the outer regions of protoplanetary disks in star forming regions. In particular they prove to be useful to identify transitional disk candidates. The goals of this work are to complement the detections of disks and the identification of transitional disk candidates in the Lupus clouds with data from the Herschel Gould Belt Survey. We extracted photometry at 70, 100, 160, 250, 350 and 500 $\mu$m of all spectroscopically confirmed Class II members previously identified in the Lupus regions and analyzed their updated spectral energy distributions. We have detected 34 young disks in Lupus in at least one Herschel band, from an initial sample of 123 known members in the observed fields. Using the criteria defined in Ribas et al. (2013) we have identified five transitional disk candidates in the region. Three of them are new to the literature. Their PACS-70 $\mu$m fluxes are systematically higher than those of normal T Tauri stars in the same associations, as already found in T Cha and in the transitional disks in the Chamaeleon molecular cloud. Herschel efficiently complements mid-infrared surveys for identifying transitional disk candidates and confirms that these objects seem to have substantially different outer disks than the T Tauri stars in the same molecular clouds.Comment: Accepted for publication in A&A. 16 pages, 9 figures, 7 table

Embedded Star Formation in the Eagle Nebula with Spitzer/GLIMPSE

We present new Spitzer photometry of the Eagle Nebula (M16, containing the optical cluster NGC 6611) combined with near-infrared photometry from 2MASS. We use dust radiative transfer models, mid-infrared and near-infrared color-color analysis, and mid-infrared spectral indices to analyze point source spectral energy distributions, select candidate young stellar objects (YSOs), and constrain their mass and evolutionary state. Comparison of the different protostellar selection methods shows that mid-infrared methods are consistent, but as has been known for some time, near-infrared-only analysis misses some young objects. We reveal more than 400 protostellar candidates, including one massive young stellar object (YSO) that has not been previously highlighted. The YSO distribution supports a picture of distributed low-level star formation, with no strong evidence of triggered star formation in the ``pillars''. We confirm the youth of NGC 6611 by a large fraction of infrared-excess sources, and reveal a younger cluster of YSOs in the nearby molecular cloud. Analysis of the YSO clustering properties shows a possible imprint of the molecular cloud's Jeans length. Multiwavelength mid-IR imaging thus allows us to analyze the protostellar population, to measure the dust temperature and column density, and to relate these in a consistent picture of star formation in M16.Comment: 16p preprint - ApJ accepte

Identification of transitional disks in Chamaeleon with Herschel

Transitional disks are circumstellar disks with inner holes that in some cases are produced by planets and/or substellar companions in these systems. For this reason, these disks are extremely important for the study of planetary system formation. The Herschel Space Observatory provides an unique opportunity for studying the outer regions of protoplanetary disks. In this work we update previous knowledge on the transitional disks in the Chamaeleon I and II regions with data from the Herschel Gould Belt Survey. We propose a new method for transitional disk classification based on the WISE 12 micron-PACS 70 micron color, together with inspection of the Herschel images. We applied this method to the population of Class II sources in the Chamaeleon region and studied the spectral energy distributions of the transitional disks in the sample. We also built the median spectral energy distribution of Class II objects in these regions for comparison with transitional disks. The proposed method allows a clear separation of the known transitional disks from the Class II sources. We find 6 transitional disks, all previously known, and identify 5 objects previously thought to be transitional as possibly non-transitional. We find higher fluxes at the PACS wavelengths in the sample of transitional disks than those of Class II objects. We show the Herschel 70 micron band to be an efficient tool for transitional disk identification. The sensitivity and spatial resolution of Herschel reveals a significant contamination level among the previously identified transitional disk candidates for the two regions, which calls for a revision of previous samples of transitional disks in other regions. The systematic excess found at the PACS bands could be a result of the mechanism that produces the transitional phase, or an indication of different evolutionary paths for transitional disks and Class II sources.Comment: Accepted for publication in A&A: 11 March 2013 11 pages, 15 figure

The dusty environment of HD 97300 as seen by Herschel and Spitzer

Aims. We analyze the surroundings of HD 97300, one of two intermediate-mass stars in the Chamaeleon I star-forming region. The star is known to be surrounded by a conspicuous ring of polycyclic aromatic hydrocarbons (PAHs). Methods. We present infrared images taken with Herschel and Spitzer using 11 different broad-band filters between 3.6 um and 500 um. We compare the morphology of the emission using cuts along different position angles. We construct spectral energy distributions, which we compare to different dust models, and calculate dust temperatures. We also derive opacity maps and analyze the density structure of the environment of HD 97300. Results. We find that HD 97300 has no infrared excess at or below 24 um, confirming its zero-age main-sequence nature. The morphology of the ring is very similar between 3.6 um and 24 um. The emission at these wavelengths is dominated by either PAH features or PAH continuum. At longer wavelengths, only the northwestern part of the ring is visible. A fit to the 100-500 um observations suggests that the emission is due to relatively warm (~26 K) dust. The temperature gradually decreases with increasing distance from the ring. We find a general decrease in the density from north to south, and an approximate 10% density increase in the northeastern part of the ring. Conclusions. Our results are consistent with the theory that the ring around HD 97300 is essentially a bubble blown into the surrounding interstellar matter and heated by the star.Comment: 6 pages, 3 figures, accepted for publication in A&

<i>Herschel</i> observations of the debris disc around HIP 92043

Context. Typical debris discs are composed of particles ranging from several micron sized dust grains to km sized asteroidal bodies, and their infrared emission peaks at wavelengths 60-100 μm. Recent Herschel DUNES observations have identified several debris discs around nearby Sun-like stars (F, G and K spectral type) with significant excess emission only at 160 μm. Aims. We observed HIP 92043 (110 Her, HD 173667) at far-infrared and sub-millimetre wavelengths with Herschel PACS and SPIRE. Identification of the presence of excess emission from HIP 92043 and the origin and physical properties of any excess was undertaken through analysis of its spectral energy distribution (SED) and the PACS images. Methods. The PACS and SPIRE images were produced using the HIPE photProject map maker routine. Fluxes were measured using aperture photometry. A stellar photosphere model was scaled to optical and near infrared photometry and subtracted from the far-infared and sub-mm fluxes to determine the presence of excess emission. Source radial profiles were fitted using a 2D Gaussian and compared to a PSF model based on Herschel observations of α Boo to check for extended emission. Results. Clear excess emission from HIP 92043 was observed at 70 and 100 μm. Marginal excess was observed at 160 and 250 μm. Analysis of the images reveals that the source is extended at 160 μm. A fit to the source SED is inconsistent with a photosphere and single temperature black body. Conclusions. The excess emission from HIP 92043 is consistent with the presence of an unresolved circumstellar debris disc at 70 and 100 μm, with low probability of background contamination. The extended 160 μm emission may be interpreted as an additional cold component to the debris disc or as the result of background contamination along the line of sight. The nature of the 160 μm excess cannot be determined absolutely from the available data, but we favour a debris disc interpretation, drawing parallels with previously identified cold disc sources in the DUNES sample

A VLA Search for Water Masers in Six HII Regions: Tracers of Triggered Low-Mass Star Formation

We present a search for water maser emission at 22 GHz associated with young low-mass protostars in six HII regions -- M16, M20, NGC 2264, NGC 6357, S125, and S140. The survey was conducted with the NRAO Very Large Array from 2000 to 2002. For several of these HII regions, ours are the first high-resolution observations of water masers. We detected 16 water masers: eight in M16, four in M20, three in S140, and one in NGC 2264. All but one of these were previously undetected. No maser emission was detected from NGC 6357 or S125. There are two principle results to our study. (1) The distribution of water masers in M16 and M20 does not appear to be random but instead is concentrated in a layer of compressed gas within a few tenths of a parsec of the ionization front. (2) Significantly fewer masers are seen in the observed fields than expected based on other indications of ongoing star formation, indicating that the maser-exciting lifetime of protostars is much shorter in HII regions than in regions of isolated star formation. Both of these results confirm predictions of a scenario in which star formation is first triggered by shocks driven in advance of ionization fronts, and then truncated approximately 10^5 years later when the region is overrun by the ionization front.Comment: 30 pages, 20 figures, 3 tables. Accepted for publication by ApJ. Full resolution figures and PS and PDF versions with full-res figures available at http://eagle.la.asu.edu/healy/preprints/hhc0

Incidence of debris discs around FGK stars in the solar neighbourhood

Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system's counterparts are the asteroid and Edgeworth-Kuiper belts. The aim of this paper is to provide robust numbers for the incidence of debris discs around FGK stars in the solar neighbourhood. The full sample of 177 FGK stars with d<20 pc proposed for the DUNES survey is presented. Herschel/PACS observations at 100 and 160 micron complemented with data at 70 micron, and at 250, 350 and 500 micron SPIRE photometry, were obtained. The 123 objects observed by the DUNES collaboration were presented in a previous paper. The remaining 54 stars, shared with the DEBRIS consortium and observed by them, and the combined full sample are studied in this paper. The incidence of debris discs per spectral type is analysed and put into context together with other parameters of the sample, like metallicity, rotation and activity, and age. The subsample of 105 stars with d<15 pc containing 23 F, 33 G and 49 K stars, is complete for F stars, almost complete for G stars and contains a substantial number of K stars to draw solid conclusions on objects of this spectral type. The incidence rates of debris discs per spectral type 0.26 (6 objects with excesses out of 23 F stars), 0.21 (7 out of 33 G stars) and 0.20 (10 out of 49 K stars), the fraction for all three spectral types together being 0.22 (23 out of 105 stars). Uncertainties corresponding to a 95% confidence level are given in the text for all these numbers. The medians of the upper limits of L_dust/L_* for each spectral type are 7.8E-7 (F), 1.4E-6 (G) and 2.2E-6 (K); the lowest values being around 4.0E-7. The incidence of debris discs is similar for active (young) and inactive (old) stars. The fractional luminosity tends to drop with increasing age, as expected from collisional erosion of the debris belts.Comment: 31 pages, 15 figures, 10 tables, 2 appendice

Dissecting the Cosmic Infrared Background with 3D Instruments

The cosmic infrared background (CIB) consists of emission from distant, dusty, star-forming galaxies. Energetically, the CIB is very important as it contains as much energy as the extragalactic optical background. The nature and evolutionary status of the objects making up the background are, however, unclear. The CIB peaks at ~150 microns, and as such is most effectively studied from space. The limited apertures of space-borne telescopes set the angular resolution that can be attained, and so even Herschel, with its 3.5m diameter, will be confusion-limited at this wavelengths at ~5mJy. The bulk of the galaxies contributing to the CIB are fainter than this, so it is difficult to study them without interferometry. Here we present the results of a preliminary study of an alternative way of probing fainter than the continuum confusion limit using far-IR imaging spectroscopy. An instrument capable of such observations is being planned for SPICA - a proposed Japanese mission with an aperture equivalent to that of Herschel and more than 2 orders of magnitude more sensitive. We investigate the potential of imaging spectrometers to break the continuum confusion limit. We have simulated the capabilities of a spectrometer with modest field of view (2'x2'), moderate spectral resolution (R~1-2000) and high sensitivity. We find that such an instrument is capable of not only detecting line emission from sources with continuum fluxes substantially below the confusion limit, but also of determining their redshifts and, where multiple lines are detected, some emission line diagnostics. We conclude that 3-D imaging spectrometers on cooled far-IR space telescopes will be powerful new tools for extragalactic far-IR astronomy.Comment: Accepted for publication in Astronomy & Astrophysic

A peculiar class of debris disks from Herschel/DUNES - A steep fall off in the far infrared

Aims. We present photometric data of debris disks around HIP 103389 (HD 199260), HIP 107350 (HN Peg, HD206860), and HIP 114948 (HD 219482), obtained in the context of our Herschel Open Time Key Program DUNES (DUst around NEarby Stars). Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 3 sigma sensitivity of a few mJy at 100 um and 160 um. In addition, we obtained Herschel/PACS photometric data at 70 um for HIP 103389. Two different approaches are applied to reduce the Herschel data to investigate the impact of data reduction on the photometry. We fit analytical models to the available spectral energy distribution (SED) data. Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths > 70 um. We investigate the significance of the peculiar shape of these SEDs and the impact on models of the disks provided it is real. Our modeling reveals that such a steep decrease of the SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a standard equilibrium collisional cascade. In contrast, a very distinct range of grain sizes is implied to dominate the thermal emission of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete and that the significance of our results depends on the version of the data reduction pipeline used. Conclusions. A new mechanism to produce the dust in the presented debris disks, deviations from the conditions required for a standard equilibrium collisional cascade (grain size exponent of -3.5), and/or significantly different dust properties would be necessary to explain the potentially steep SED shape of the three debris disks presented. (abridged)Comment: 14 pages, 4 figures, accepted by A&